Measurement of the Generalized Polarizabilities of the Proton in Virtual Compton 
Scattering at (5^=0.92 and 1.76 GeV^ : I. Low Energy Expansion Analysis 



0^ 



S. Jaminion,^ N. Degrande,^ G. Laveissiere,^ C. Jutier,^'^ L. Todor,'^ R. Di Salvo, ^ L. Van Hoorebeke,^ L.C. Alexa,^ 
B.D. Anderson,^ K.A. Amol,^ K. ArundellJ G. Audit,^ L. Auerbach,^ F.T. Baker/" M. Baylac,^ J. Berthot/ 
P.Y. Bertin/ W. Bertozzi," L. Bimbot/^ W.U. Boeglin,!^ E.J. Brash,-* V. Breton,^ H. Breuer,*-* E. Burtin,^ 
J.R. Calarco,i5 L.S. Cardman,!^ C. Cavata,^ C.-C. Chang," J.-P. Chen,!^ E. Chudakov,!^ E. Cisbani,!^ 

D.S. Dale,!^ C.W. de Jager,!*^ R. De Leo,*^ A. Deur,^-^^ N. d'Hose,^ G.E. Dodge,^ J.J. Domingo,!'^ 
L. Elouadrhiri,!^ M.B. Epstein,*^ L.A. Ewell,*'' J.M. Finn,"^ K.G. Fissum,^! H. Fonvieille,i G. Fournier,^ 
B. Frois,8 S. FruUani,!^ C. Furget,20 H. Gao,"'^! j. Gao," F. Garibaldi," A. Gasparian,22. is S. Gilad," 
R. Gilman,23. 16 A. Glamazdin,24 C. Glashausser,23 J. Gomez,!^ V. Gorbcnko,^^ p. Grenier,* P.A.M. Guichon,^ 
J.O. Hansen,i6 R. Holmes,^^ M. Holtrop,!^ C. Howell,2i G.M. Huber,^ C.E. Hyde-Wright,^ S. Incerti,^ M. lodice," 
J. Jardillier,8 M.K. Jones,^-!^ W. Kahl,^^ S. Kato,^^ A.T. Katramatou,^ J.J. Kelly,!"* Kerhoas,** A. Ketikyan,^^ 
^ . M. Khayat,^ K. Kino.^^ S. Kox,20 L.H. Kramer.^^ K.S. Kumar,^^ G. Kumbartzki,23 M. Kuss,i^ A. Leone,^" 

O ; J- J- LeRose,**^ M. Liang,*^ R.A. Lindgren,3i N. Liyanage,"'^! qj Lolos,^ R.W. Lourie,^^ R. Madey,^ K. Maeda,^^ 
CN . S. Malov,23 D.M. Manley,^ G. Marchand,*^ D. Marchand,^ D.J. Margaziotis,*^ P. Markowitz,*^ J. Marroncle,^ 

C ; J. Martino,^ C.J. Martoff,^ K. McCormick,3'23 j. Mclntyre,^^ S. Mehrabyan,^^ F. Merchez,20 Z.E. Meziani,^ 

^ ■ R. Michaels,!^ C.W. Miller,^^ J.Y. Mougey,20 S.K. Nanda,!^ D. Neyret,^ E.A.J.M. Offermann,!^^ Z. Papandreou,"^ 

i G.F. Perdrisat,"^ R. Perrinc^o G.G. Petratos,^ S. Platchkov,^ R. Pomatsalyuk,^'! D.L. Prout,^ V.A. Punjabi,33 

^ T. Pussieux,^ G. Quemener,'^'^^ R.D. Ransome,^^ Q. Ravel,* J.S. Real,^" F. Renard,^ Y. Roblin,*'*^ D. Rowntree," 

G. Rutledge,^ P.M. Rutt,^^ A. Saha,*^ T. Saito,^^ A.J. Sarty,34 A. Serdarevic,*- *6 T. Smith,*^ G. Smirnov,* 
CN ; K. Soldi,35 P. Sorokin,24 P.A. Souder,^^ R. Suleiman,^'" J.A. Templon,*" T. Terasawa,^^ R. Tieulent,^" 

1> '. E. Tomasi-Gustaffson,^ H. Tsubota,^^ H. Uenc^e P.E. Ulmer,^ G.M. Urciuoh,*^ R. Van De Vyver,^ R.L.J. Van dcr 
j^ggj. 4,16 Vernin,^ B. Vlahovic,^^ H. Voskanyan,^^ E. Voutier,^" J.W. Watson, ^ L.B. Weinstein,^ 
■ K. Wijesooriya,'^ R. Wilson,^^ B.B. Wojtsekhowski,**^ D.G. Zainea,^ W-M. Zhang,^ J. Zhao,** and Z.-L. Zhou** 
<N ! (The Jefferson Lab Hall A Collaboration) 

I 'Universite Blaise Pascal/IN2P3, F-63177 Aubiere, France 

. ^ University of Gent, B-9000 Gent, Belgium 

--^ ■ ^Old Dominion University, Norfolk, VA 23529 

\ ^University of Regina, Regina, SK S4S 0A2, Canada 

^Kent State University, Kent OH U242 
^ ' ''California State University, Los Angeles, CA 90032 

I ''College of William and Mary, Williamsburg, VA 23187 

^ ! '^CEA Saclay, F-91191 Gif-sur-Yvette, France 

• Temple University, Philadelphia, PA 19122 

.5^ \ "'University of Georgia, Athens, GA 30602 

. '^'^ Massachusetts Institute of Technology, Cambridge, MA 02139 

' ^^Institut de Physique Nucleaire, F-91406 Orsay, France 

Florida International University, Miami, FL 33199 
University of Maryland, College Park, MD 20742 
University of New Hampshire, Durham, NH 03824 
Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 
'''INFN, Sezione Sanita and Istituto Superiore di Sanitd, 00161 Rome, Italy 
University of Kentucky, Lexington, KY 40506 
^^INFN, Sezione di Bari and University of Bari, 70126 Bari, Italy 
^" Laboratoire de Physique Subatomique et de Cosmologie, F-38026 Grenoble, France 
^^Duke University, Durham, NC 27706 
Hampton University, Hampton, VA 23668 
^'^ Rutgers, The State University of New Jersey, Piscataway, NJ 08855 
Kharkov Institute of Physics and Technology, Kharkov 61108, Ukraine 
Syracuse University, Syracuse, NY 13244 
^^Yamagata University, Yamagata 990, Japan 
^''Yerevan Physics Institute, Yerevan 375036, Armenia 
^^Tohoku University, Sendai 980, Japan 
^^Princeton University, Princeton, NJ 08544 
^"INFN, Sezione di Lecce, 73100 Lecce, Italy 
^' University of Virginia, Charlottesville, VA 22901 
'^^ State University of New York at Stony Brook, Stony Brook, NY 11794 



2 



Norfolk State University, Norfolk, VA 23504 
Florida State University, Tallahassee, FL 32306 
North Carolina Central University, Durham, NC 27707 
•'^Harvard University, Cambridge, MA 02138 

Virtual Compton Scattering is studied at the Thomas Jefferson National Accelerator Facility at 
low Center-of-Mass energies, below pion threshold. Following the Low Energy Theorem for the 
ep ep7 process, we obtain values for the two structure functions Pll — Ptt /e and Plt at 
four-momentum transfer squared Q'^= 0.92 and 1.76 GeV"^. 

PACS numbers; 13.60.-r,13.60.Fz 



The electric and magnetic polarizabilities of the nu- 
cleon reflect its response to a static electromagnetic field. 
These are fundamental observables of the ground state, 
closely related to the entire excitation spectrum of the nu- 
cleon. The polarizabilities of the proton have been mea- 
sured in Real Compton Scattering (RCS) experiments 
7P — > 7p; see e.g. Contrary to atomic polarizabilities, 
which are of the size of the atomic volume , the proton 
electric polarizability is much smaller than one cubic fm, 
the volume scale of a nucleon. In a simplified harmonic 
oscillator model, such a small electric polarizability is a 
natural indication of the intrinsic relativistic character of 
the nucleon. The smallness of the proton magnetic polar- 
izability (3m relative to a^; reflects a strong cancellation 
of para- and dia-magnetism in the proton. 

A theoretical study of the Virtual Compton Scattering 
(VCS) reaction 7*p — > 7p, at threshold but at arbitrary 
virtuality of the initial photon, led to the concept of 
Generalized Polarizabilities (CP) The CPs measure 
the spatial variation of the polarization of the proton, in- 
duced by external electric- or magnetic-dipole fields as a 
function of Q^. They provide an original way to probe 
the nucleon structure. After the NE-18 experiment 3 
and the pioneering VCS experiment at MAMI the 
E93-050 experiment Q was performed at higher energies 
at Jefferson Lab ( JLab) . The studied reaction channel is 
the exclusive photon clcctroproduction ep — > ep^. In this 
Letter, we report the results of the low-energy expansion 
(LEX) analysis of this experiment, based on the devia- 
tions of the cross section with respect to the Low Energy 
Theorem. In a companion Letter (referred to as II) we 
report the results of a Dispersion Relation analysis of our 
data 0. 

According to P. Guichon et al. the unpolarized 
cross section for the reaction ep — > ep7 at small mo- 
mentum q' of the final photon can be written: 

^-0 = vi{Pll--Ptt) + viPlt (1) 
e 

where </>, vi, V2 are kinematical coefficients, q' is the mod- 
ulus of the three-momentum of the final photon in the 
(7p) CM frame and e is the virtual photon polarization. 
The differential elements of dPa are the scattered electron 
momentum and solid angle in the lab frame and the pro- 
ton solid angle in the (jp) CM frame. vI/q contains the 
effect of the CPs through the structure functions Pll, 



Ptt, Plt- ^'o represents the effect of CPs to first or- 
der in q' , or equivalently the truncation to electric- or 
magnetic-dipole radiation. It has been shown to contain 
six independent (dipole) CPs 0, , combined into the 
structure functions. (f'crBH+Born corresponds to the co- 
herent sum of the Bethe-Heitler (BH) and the VCS Born 
amplitudes. It depends only on the elastic form factors 
of the proton and is a particular case of Low's low-energy 
theorem [Toj for threshold photon production. 

The goal of the measurement, performed at fixed e, is 
to extract the two structure functions Pll — Ptt 1^ and 
Plt by extrapolating the quantity 

(dV^^^ - d5a^^+^°™)/(gV) (2) 

to g' = 0. This is done at a fixed value of g, the CM 
three-momentum of the VCS virtual photon. Thus, for 
the small but non- vanishing values of g' that we consider, 
the corresponding squared momentum = 9^ ~ 9o 
pends on g' (through the CM virtual photon energy go)- 
Though this dependence is very weak for our kinemat- 
ical conditions, we point out that in the following, 
will refer to the value corresponding to g' = 0, that is 

Q2 ^ 2A/p • {^Ml +q'^-Mp). Without hsting the com- 
plete expression of the structure functions in terms of the 
CPs la, Ell, we just mention that Pll is proportional to 
the electric CP asiQ^) and that Plt contains the mag- 
netic GP (3m {Q"^)- When tends to zero, these two 
CPs give the well-known electric and magnetic polariz- 
abilities obtained in RCS. The structure function Ptt is 
a combination of spin- flip GPs. 

The apparatus and running conditions of the JLab ex- 
periment are detailed elsewhere [lllTllTl IlllTlliTt . 
Here a brief summary is given. An electron beam of 
4.030 GeV energy was directed onto a 15 cm liquid hy- 
drogen target. The two Hall A Spectrometers were used 
to detect the scattered electron and the outgoing proton 
in coincidence, allowing the identiflcation of the exclu- 
sive reaction ep epj by the missing-mass technique. 
This experiment makes use of the full capabilities of the 
accelerator and the Hall A instrumentation Ts]: 100% 
duty cycle, high resolution spectrometers, high luminosi- 
ties (up to 4 X 10^^ cm^^ • s^^). 

The data for the LEX analysis are divided into two sub- 
sets corresponding to data taking in two different ranges 
of [0.85, 1.15] and [1.6, 2.1] GeV^. Spectrometer 
optics are optimized for the experiment, and dedicated 
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procedures allow adjustment of the main experimental 
offsets in energies, angles and positions. For each event, 
variables such as g', or the CM polar and azimuthal an- 
gles 9 and (p of the outgoing photon w.r.t. (f, are obtained 
by reconstructing the missing particle. The acceptance 
calculation is provided by a dedicated Monte-Carlo simu- 
lation 19] including resolution effects. Radiative correc- 
tions are applied following the exponentiation method 
of Ref . , the acceptance-dependent part being imple- 
mented in the simulation. The acceptance calculation 
requires a realistic shape for the sampling cross section. 
For this reason, simulated events are generated using the 
(BH-|-Born) cross section plus a first-order CP effect, in- 
troduced iteratively. Clean event samples are obtained by 
applying the event selection method. It involves simple 
cuts, e.g. a window around in missing mass squared to 
select the photon peak, or the condition q' < 126 MeV/c 
to stay below pion threshold. Other cuts are less triv- 
ial; special attention was paid to obtain a well-defined 
acceptance and to eliminate protons punching through 
the spectrometer entrance collimator. The same cuts 
are applied to experimental and simulated events. The 
absolute normalization of the cross section is obtained 
from the knowledge of the beam charge, target density, 
and detectors efficiencies. More details can be found in 
Refs. [Illlllll. 

For each of the two data sets, the photon electropro- 
duction cross section is determined at a fixed value of q 
(1.080 and 1.625 GeV/c) and e (0.95 and 0.88). Events 
are binned in q' (30 MeV/c wide), and in 9 and ip (see 
Fig.Q. The advantage of the JLab experiment is that it 
produces a large out-of-plane acceptance, thanks to the 
Lorcntz boost (CM — > lab) focusing the outgoing proton 
around q. An example of out-of-plane cross-section data 
is given in Fig.[T] At the lowest q' of 45 MeV/c, the CP 
effect is very small (< 3% of dPa) and not visible on the 
plot; at the highest q' below pion threshold (105 MeV/c) 
it reaches 10-15% . We confirmed the central value of the 
absolute normalization of the experiment to the 1% level 
in the lowest g'-bin, where the experimental cross sec- 
tion is close to the known (BH-|-Born) cross section. For 
this test and throughout our analyses, the (BH-|-Born) 
cross section is calculated using the parametrization of 
Ref. ,2]J for the proton EM form factors. 

The method to extract the structure functions is de- 
duced from Eq. ^ For each bin in {9, if), one measures 
^b^EXP g^-j- ggveral finite values of q' , and extrapolates 
the quantity AM = {d^a^^^ - dV^^+-»°™)/(09') to 
q' = 0, yielding the value of ^'q. In our data, does 
not exhibit any significant g'-dependence, so the extrap- 
olation to g' = is done in each bin in (0, ip) by averaging 
AA^ over q' . The resulting ^0 term is then fitted as a 
linear combination of two free parameters, which are the 
structure functions Pll ~ PttI^ and P^t (see Fig. 
Their statistical error is given by the fit, and their system- 
atic error is calculated from four sources added quadrat- 
ically: 1) ± 2 MeV on beam energy, 2) ± 0.5 mrad on 
horizontal angles, 3) ± 2.3% on overall absolute cross 
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FIG. 1: {ep —* epy) cross section for the lowest and highest 
q' bin, at 40° out-of-plane (including symmetry w.r.t. the 
lepton plane). Only statistical errors are shown. Here the 
out-of-plane angle (or latitude) is the polar angle of the out- 
going photon when the polar axis is chosen perpendicular to 
the lepton plane. The abscissa is the azimuthal angle (or lon- 
gitude) using this convention; negative values are defined as 
in Ref. |5|. The full curve is the (BH-fBorn) cross section, 
the dashed curve includes the first-order GP effect (Eq.0 as 
fitted in this analysis. 



section normalization and 4) ± 2% due to possible cross 
section shape distortions. The value of the reduced of 
the fit (1.22 and 1.50, cf. Fig. is one way to express 
that the LEX holds reasonably well at our kinematics. 

Our results, presented in Table together with the 
previous results at lower Q^, show a strong fall-off of the 
structure functions with momentum transfer, similarly to 
form factors. This behavior is expected since CPs can be 
seen as "form factors of a nucleon embedded in an EM 
field". From a comparison with theoretical predictions, 
one should point out that most model calculations of the 
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FIG. 2: A graphical representation of the LEX fit (straight 
line) for each data set. Black circles correspond to out-of- 
plane data, and the inner plot is a zoom on the lepton plane 
data (triangles). ^o,vi and V2 are defined in the text. 
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TABLE I: Compilation of the VCS structure functions. The 
first hne corresponds to the RCS result 0. In all cases the 
first error is statistical, and the second one is the total sys- 
tematic error. 



Ref. 




e Pll — Ptt/i^ 


Plt 




(GeV^) 


(GeV-2) 


(GeV-2) 


Previous experiments 


UJ 





81.3 


-5.4 






± 2.0 ± 3.4 


± 1.3 ± 1.9 




0.33 


0.62 23.7 


-5.0 






± 2.2 ± 4.3 


± 0.8 ± 1.8 


This experiment 




0.92 


0.95 1.77 


-0.56 






± 0.24 ± 0.70 


± 0.12 ± 0.17 




1.76 


0.88 0.54 


-0.042 






± 0.09 ± 0.20 


± 0.047 ± 0.055 



GPs, e.g. Refs. HHHHIl are valid at much smaller 
than the ones covered in this experiment. One notable 
exception is the Dispersion Relation (DR) model [25^ : 
our extraction of VCS structure functions in the same 
experiment, based on the DR approach, is the subject of 
a companion Letter ■ The results presented here should 
stimulate theoretical calculations of GPs at intermediate 
and high Q"^ . It would be of great interest if the models 



that predict the form factors Ge and Gm at > 1 
GeV^ could give predictions of the GPs in the same 
range. 

In summary, we have studied the process ep —>■ epj 
using the unique capabilities of the JLab electron beam 
and the Hall A instrumentation. We have measured the 
two unpolarized structure functions of Virtual Compton 
Scattering Pll — Ptt/^ and Plt, providing a new piece 
of information on the electromagnetic structure of the 
nucleon, namely on the behavior of the proton General- 
ized Polarizabilities in the range ~ 1-2 GeV^. Other 
VCS experiments at low energy will measure these struc- 
ture functions (MIT-Bates [26|, MAMI ^) and double 
polarization is foreseen j23| as a way to separate the six 
lowest order GPs. 
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